Phase Transition Of The Spin-One Square-Lattice Anisotropic Heisenberg Model 
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The coupled cluster method (CCM) is applied to the spin-one anisotropic Heisenberg antiferro- 
magnet (HAF) on the square lattice at zero temperature using a new high-order CCM ground-state 
formalism for general quantum spin number (s > 1/2). The results presented constitute the first 
such application of this new formalism, and they are shown to be among the most accurate re- 
sults for the ground-state energy and sublattice magnetisation of this model as yet determined. We 
"track" the solution to the CCM equations at a given level of approximation with respect to the 
anisotropy parameter, A, from the trivial Ising limit (A — > oo) down to a critical value A c , at which 
point this solution terminates. This behaviour is associated with a phase transition in the system, 
and hence a primary result of these high-order CCM calculations is that they provide an accurate 
and unbiased (i.e., ab initio) estimation of the position of the quantum phase transition point as a 
function of the anisotropy parameter. Our result is, namely, that this point occurs at (or slightly 
below) the isotropic Heisenberg point at A = 1, for this model. 

PACS numbers: 75.10.Jm, 75.50.Ee, 03.65.Ca 
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The area of quantum magnetic insulating systems at 
zero temperature has become increasingly well under- 
stood for one-dimensional (or quasi-one-dimensional) lat- 
tices via the existence of well-known exact solutions such 
as the Bethe ansatz (BA) |l], |[ and also via more re- 
cent density matrix renormalisation group (DMRG) cal- 
culations H ^, ^). Similarly, zero-temperature quan- 
tum Monte Carlo (QMC) calculations jg, [jj have been 
shown to yield very accurate results for spin-half, two- 
dimensional systems. QMC techniques, however, are 
severely limited in their range of application by the pres- 
ence of the infamous "sign problem," although we note 
that for non-frustrated spin systems it is often possi- 
ble to determine a "sign rule" ^ which completely 
circumvents this minus-sign problem. By contrast, the 
technique of quantum many-body theory known as the 
coupled cluster method (CCM) is neither limited by the 
presence of frustration nor by_ the spatial dimensional- 
ity of the lattice 0, ^ |f, p 



. Indeed, a great 

strength of this method is that it is able to determine with 
great accuracy the positions of quantum phase transition 
points of quantum systems within an ab initio, and thus 
essentially unbiased, framework. These quantum phase 
transitions typically arise as some parameter within the 
Hamiltonian is varied, thus driving the system from one 
phase into another. 

The behaviour of the quantum Heisenberg antiferro- 
magnet (HAF) on the linear chain is highly dependent 
on whether the quantum spin number, s, has an integer 
or a half-odd- integer value. Indeed, Haldane 16 first 
predicted that the spin-one isotropic HAF would contain 
an excitation energy gap, and this prediction has subse- 
quently and conclusively been shown numerically to be 
correct by exact diagonalisations of short chains |17]] and 
by more recent and extremely accurate DMRG calcu- 
lations ||. We note that this is in stark contrast to the 
exact BA solution p], || of the spin-half isotropic HAF for 
the linear chain, which contains no such gap. However, 



such exact diagonalisations and DMRG calculations are 
considerably more difficult for systems of higher spatial 
dimensionality. In this article, we focus on the specific 
case of the spin-one anisotropic HAF on the square lattice 
at zero temperature using the CCM, especially in relation 
to its quantum phase transition, by applying a high-order 
CCM formalism for general quantum spin number for the 
first time. 

The spin-one anisotropic HAF is given by, 



H = 



(1) 



where the symbol indicates nearest-neighbour 

bonds on the square lattice, and where each bond is 
counted once and once only. For A > 1, the ground state 
of this model contains non-zero Neel-like order. The pre- 
cise value of the phase transition point, at which this 
ordering breaks down is not exactly known although it 
is believed to be at (or near to) the Heisenberg point 
at A = 1. Very accurate results for the values of various 
ground-state properties of this model have been obtained 
via spin-wave theory (SWT) and cumulant series ex- 
pansions |lj| . Both sets of results indicate that approxi- 
mately 80% of the classical ordering remains for the spin- 
one, square-lattice HAF model at A = 1. We note how- 
ever that cumulant series expansions make an explicit as- 
sumption that the position of the phase transition point is 
exactly at A = 1 in order to perform Pade resummations 
of the otherwise divergent perturbation series. Also, con- 
ventional SWT essentially appears to "build in" a phase 
transition point at A = 1 for any lattice, at which point 
the excitation spectrum becomes soft. This is, of course, 
now known to be incorrect J|, [l6|, 17 for the spin-one 
linear-chain HAF model, as noted above. 

We now briefly describe the general CCM formalism, 
and the interested reader is referred to Refs. [[lO], [ll], 
HI pi pi 111 for further details. The exact ket and 
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FIG. 1: CCM results for the ground-state energy per spin of 
the spin-one anisotropic Heisenberg antiferromagnet on the 
square lattice using the LSUBm and SUBm-m approxima- 
tions with m = {2, 4, 6}. The boxes indicate the CCM critical 
points. 



FIG. 2: CCM results for the sublattice magnetisation of 
the spin-one anisotropic Heisenberg antiferromagnet on the 
square lattice using the LSUBm and SUBm-m approxima- 
tions with m = {2, 4, 6}. 



bra ground-state energy eigenvectors, |\&) and (^l, of a 
many-body system described by a Hamiltonian H are 
parametrised within the single-reference CCM as follows: 

|v&)=e s |<i>) ; S = Y,SiCf: 

IjtO 

(*| = ($|£e- s ; 5=1 + ^5/(77- (2) 

The single model or reference state <&} is required to have 
the property of being a cyclic vector with respect to two 
well-defined Abelian subalgebras of multi-configurational 
creation operators {Cj} and their Hcrmitian-adjoint de- 
struction counterparts {Cj = (Cj)' 1 '}. Thus, |$) plays 
the role of a vacuum state with respect to a suitable set 
of (mutually commuting) many-body creation operators 
{Cj}. Note that Cj\§) = 0, V I ^ 0, and that Cq = 1, 
the identity operator. 

We note that the exponentiated form of the ground- 
state CCM parametrisation of Eq. (||) ensures the correct 
counting of the independent and excited correlated many- 
body clusters with respect to |$) which are present in 
the exact ground state |^). It also ensures the exact 
incorporation of the Goldstone linked-cluster theorem, 
which itself guarantees the size-extensivity of all relevant 
extensive physical quantities. 

The CCM equations are defined by the following cou- 
pled set of equations, 

<i>|C7e- s #e s |$) = 0, V 1^0 ; (3) 
(1>\Se- s [H,C+}e s \<S>) = 0, V I ^ . (4) 



We furthermore note that the ground-state energy at the 
stationary point has the simple form 

E g = E g ({S!}) = {<S>\ C - s Hc s \<S>) . (5) 

It is important to realize that this (bi-)variational formu- 
lation does not lead to an upper bound for E g when the 
summations for S and S in Eq. (^|) are truncated, due 
to the lack of exact Hermiticity when such approxima- 
tions are made. However, it is clear that the important 
Hellmann-Feynman theorem is still preserved in all such 
approximations. 

The CCM formalism is exact in the limit of inclusion 
of all possible multi-spin cluster correlations for S and 
S, although in any real application this is usually impos- 
sible to achieve. It is therefore necessary to utilise var- 
ious approximation schemes within S and S. The three 
most commonly employed schemes previously utilised 
have been: (1) the SUBro scheme, in which all correla- 
tions involving only n or fewer spins are retained, but no 
further restriction is made concerning their spatial sepa- 
ration on the lattice; (2) the SUBn-m sub-approximation, 
in which all SUBn correlations spanning a range of no 
more than m adjacent lattice sites are retained; and (3) 
the localised LSUBm scheme, in which all multi-spin cor- 
relations over distinct locales on the lattice defined by m 
or fewer contiguous sites are retained. We note that in 
order to carry out such LSUBm or SUBn-m calculations 
to high order in the truncation indexes n and m we must 
rely on computer-algebra techniques in order to gener- 
ate the corresponding sets of coupled equations. These 
computational techniques are based on the fact that the 
similarity transformed Hamiltonian in Eqs. (&@) may 
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TABLE I: CCM results for the ground state of the spin-one 
Heisenberg antiferromagnet at A = 1 on the square lattice 
using the LSUBm approximation with m = {2,4,6}. Values 
for the CCM critical points, A c , of the anisotropic model as 
a function of the anisotropy, A, are also presented. Note that 
Nf indicates the number of fundamental clusters at each level 
of approximation. 
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TABLE II: CCM results for the ground state of the spin-one 
Heisenberg antiferromagnet at A = 1 on the square lattice 
using the SUBm-m approximation with m — {2, 4, 6}. Values 
for the CCM critical points, A c , of the anisotropic model as 
a function of the anisotropy, A, are also presented. Note that 
Nf indicates the number of fundamental clusters at each level 
of approximation. 
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a See Ref. 
6 See Ref. 
c See Ref. 
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be written purely in terms of creation operators acting 
on the model state, so that evaluation of Eq. (^) be- 
comes purely a matter of pattern matching. The inter- 
ested reader is referred to Refs. [^2j and for a full 
explanation of how this is achieved for spin-half models. 
We note that the generalisation of this procedure to the 
case of arbitrary quantum spin number, s, is relatively 
straightforward, although detailed, and will be described 
fully elsewhere pfjfl . 

We now apply the CCM formalism outlined above to 
the specific case of the anisotropic HAF, and we choose 
the Neel state, in which the spins lie along the z-axis, 
to be the model state. Furthermore, we perform a ro- 
tation of the local axes of the up-pointing spins by 180° 
about the y axis, so that spins on both sublattices may 
be treated equivalently. The (canonical) transformation 
is described by, 

s x -s x , s« -> s v : s z -> -s z . (6) 

The model state now appears mathematically to consist 
of purely down-pointing spins in these rotated local axes. 
In terms of the spin raising and lowering operators sj: = 
s'f, ± is| the Hamiltonian may be written in these local 
axes as , 

H =-- 4 T l { s t4+ s i s J+ 2AsZ ^} > ( 7 ) 

(id) 1 J 

where i and j runs over all nearest neighbours, although 
each nearest-neighbour bond is counted once and once 
only. Furthermore, the sublattice magnetisation, M, (af- 
ter rotation of the local spin axes) is given by, 

1 N 

M = -nH< ■ (8) 



In the limit A — > oo we note that the model state is the 
exact ground eigenstate of the Hamiltonian of Eq. (Q). 
Hence, all of the CCM correlation coefficients are zero at 
this point. We may track this solution for decreasing val- 
ues of A until we reach a critical point, A c , at which point 
the real solution to our CCM equations for the LSUBm 
and SUBm-m approximation schemes terminates. This 
is associated with a phase transition in the real system 
jl2|, P, [F§|n|. Similar behaviour was also seen for this 
model for the CCM SUB2 approximation for this model 
(see Ref. @), and once more this was associated with 
a phase transition of the real system. 

We note that the "raw" results for the ground-state 
energy, sublattice magnetisation, and critical points, A c , 
have been obtained using both the LSUBm and SUBm- 
m approximation schemes. For both approximation 
schemes we may extrapolate these results as a function 
of m in the limit m — ► oo in order to obtain even better 
results. We note that SUB2-m calculations and the full 
SUB2 calculation have also previously been performed 
p"j"| for the anisotropic square-lattice HAF with general 
quantum spin number. It was therefore possible to de- 
termine accurately the manner in which SUB2-m results 
scale as a function of m, namely, that as m — ► oo the 
ground-state energy and the critical points, E g (m)/N 
and A c (m) respectively, scale linearly with 1/m 2 , and 
the sublattice magnetisation, M(rn), scales linearly with 
1/m. By analogy, we utilise a similar procedure here in 
order to extrapolate raw LSUBm and SUBm-m results. 
However, as we may utilise LSUBm and SUBm-m results 
for m = {2, 4, 6} only, a quadratic function is fitted to 
these data in order to obtain the best possible results. A 
full and comprehensive explanation of the extrapolation 
process of CCM LSUBm expectation values is given in 
Ref. [jl5| for the spin- half XXZ model for a variety of lat- 
tices, and the interested reader is referred to this article 



4 



for further information. 

CCM results for the position of the critical point using 
the LSUBm and SUBm-m approximation schemes are 
presented in Tables | and |J. The extrapolated result of 
A c = 0.98 for both the LSUBm and for the SUBm-m 
approximation scheme indicates that the phase transi- 
tion point is at (or perhaps slightly below) the isotropic 
Heisenberg point (A = 1). The strength of the CCM is 
that it can provide such an accurate value for the position 
of the quantum phase transition point as a function of the 
anisotropy within an ab initio, and thus fully unbiased, 
framework. 

Figure [l] and Tables | and |n] indicate that the results 
for the ground-state energy of the spin-one HAF con- 
verge extremely quickly with m for both the LSUBm 
and SUBm-m schemes. Indeed, even the "raw" un- 
extrapolated results provide excellent estimates of the 
ground-state energy, although the extrapolated results of 
Eg/N — —2.3292 for the LSUBm approximation scheme 
and Eg/N — —2.3291 for SUBm-m approximate scheme 
are certainly even more accurate. By way of compari- 
son, we note that third-order SWT gives a result of 
Eg/N — —2.3282 and cumulant series expansions jl9| a 
result of Eg/N = -2.3279(2). 

CCM results for the sublattice magnetisation are found 
to be similarly well converged as a function of the trunca- 
tion index m for both the LSUBm and SUBm-m schemes, 
as indicated in Fig. |[ The extrapolated results of 
M = 0.8049 and M = 0.8067 for the LSUBm and SUBm- 
m schemes, respectively, are again in excellent agreement 



with the results of third-order SWT [fil| and cumulant 
series expansions [|l9|, which give values of M = 0.8043 
and M = 0.8039(4) respectively. 

In this article it has been shown that the coupled 
cluster method (CCM) provides quantitatively accurate 
results for the ground-state properties of the spin-one 
square-lattice isotropic HAF by comparison with results 
of third-order spin- wave theory (SWT) and cumulant se- 
ries expansions. Furthermore, we note that the results 
presented in this article constitute the first application 
of the CCM using a new high-order CCM ground-state 
formalism for general (s > 1/2) quantum spin num- 
ber. (A fuller account of this new high-order formalism 
will be published elsewhere [^oft .) The best estimates 
of the ground-state energy of this model were found to 
be Eg/N = —2.3292 using the LSUBm approximation 
scheme and E g /N = —2.3291 using the SUBm-m approx- 
imation scheme via heuristic extrapolation to the (exact) 
limit m — > oo. The best estimates of the sublattice mag- 
netisation are M — 0.8049 using the LSUBm approxi- 
mation scheme and M — 0.8067 using the SUBm-m ap- 
proximation scheme, again via extrapolation to the limit 
m — > oo. The most important result of the CCM calcula- 
tions presented here for this model is the prediction that 
the phase transition point of the spin-one square-lattice 
anisotropic HAF is at (or slightly below) the isotropic 
Heisenberg point. A strength of the method is that this 
prediction may be made using an unbiased (i.e., ab initio) 
treatment. 
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